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ABSTRACT 



Aims. We improve the earlier restrictions on parameters of the dark matter (DM) in the form of a sterile neutrino. 
Methods. The results were obtained from non-observing the DM decay line in the X-ray spectrum of the Milky Way (using the recent 
XMM-Newton PN blank sky data). We also present a similar constraint coming from the recent XMM-Newton observation of Ursa 
Minor - dark, X-ray quiet dwarf spheroidal galaxy. 

Results. The new Milky way data improve on (by as much as the order of magnitude at masses ~ 3.5 keV) existing constraints. 
Although the observation of Ursa Minor has relatively poor statistics, the constraints are comparable to those recently obtained 
using observations of the Large Magellanic Cloud or M31. This confirms a recent proposal that dwarf satellites of the MW are very 
interesting candidates for the DM search and dedicated studies should be made to this purpose. 

Key words. Galaxy: halo - Galaxies: individual: Ursa Minor - Cosmology: dark matter - X-rays:galaxies 
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This past year has seen a lot of activity, devoted to 
searching for the decay signals of the DM particle in X- 
ray sp ectra of various astrophysical ob j ects ( Boyarsky et alj 
l2006aflbllcl iRiemer-Sdrensen et alj 120061 IWatson et alj [2006: 
Riemer-S0rensen et al. 2006, e.g). It was noticed long ago 
bv lDodelson & Wid row ( 1994) that a right-handed neutrino with 
masses in the keV range presents a viable warm dark matter 
(WDM) candidate. Such a particle possesses a specific radia- 
tive decay channel, so one can search for its decay line in the 
X-ray spectra of astrop hysical objects ( Dol gov & Hansenll2002t 
lAbazajian et~ai1l2001bl) . 

Recently, the interest in the sterile neutrino as a DM can- 
didate has been greatly rev italized. First, the discove ry of neu- 
trino oscillations (see e.g. Strumia & Vissani (2006) for a re- 
view) strongly suggest the existence of right-handed neutrinos. 
Probably the easiest way to explain the data on oscillations is by 
adding several right-handed, or sterile, neutri nos, to the Standard 
Mod el. It has been demonstrated r ecently in lAsaka et al.l d2005l) 
and lAsaka & Sh aposhnikov ( 2005) that a simple extension of the 
Standard Model by three singlet fermions with ma sses smaller 
than the electroweak scale (dubbed the vMSM in Asak a et ail 
120051) allows accommodation of the data on neutrino masses 
and mixings, allows baryon asymmetry of the Universe to be 



explained, and provides a candidate for the dark matter particle 
in the form of the lightest of the sterile neutrinosQ. 

Secondly, warm DM with the mass of particle in keV range 
can ease the problem of the dark halo st ructures in c ompar - 
ison with the cold dark matter scenario iBode et alJ §001); 
iGoerdt et alj d2006l) . By determining the ma tter power spectrum 
fro m the Lyman-q fo rest data from SDSS ISeliak et"aU {2006) 
and IViel et alj d2006l) argue that the mass of the DM parti- 
cl es should be in th e range > 10 keV (> 14 keV in the case 
of lSeliak et aIB2 006). As this method gives direct bounds for the 
free-streaming length of the neutrinos, the bounds on the mass 
of the DM particle depend on the momentum distribution func- 
tion of the sterile neutrinos and, therefore, on their production 
mechanism. The results quoted abov e are claimed for the sim- 
plest Dodelson- Widrow model dl994l) . 

At th e same time, stud i es of the Fornax dw arf spheroidal 
galaxy dGoerdt et alj 120061: IStrigari et ail [2006) disagree with 
th e predictions of C DM mo dels and sugge st lower mass than 
in ISeliak et alj d2006l) and IViel et all (2006) for the DM parti- 
cle Mqm ~ 2 keV. This result agrees with the earlier studies 
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1 The vMS M does not explai n the unconfirmed results of the LSND 
experiment ( Aguilar et al. 2001). There are other models that try to ac- 
count for it by introducing a sterile neutrino with the mass around 1 
eV. There are also models that explain not all, but only some of these 
p henomena (e.g. LS ND and DM, but not the baryon asymmetry as e.g. 
in de Gouvea 2005) We do not give any review here. We would like to 
stress that, although our work is motivated by the recent results on the 
vMSM, our method and results do not rely on any particular model. 
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of lHansen et all d2002l) and lViel et all (120051) . which used a dif- 
ferent dataset. For other interesting applica tions of the steril e 
neutrinos with the ma s s ~ few keV see e.g . Kusenko (2006b), 
iBiermann & Kusenkol (120061). [Stasielak et all d2006). Ku senkol 
d2006al) . and lHidaka & Fuller! (120061). 



Riemer-S0rensen et al 



gued in 

U (120061) 



It has been argued in iBovarskv et al.1 d2006c) and 



that the preferred targets for 



observations are objects from the local halo, includin g our 
own M ilky Way and its satellites. In particular, Boyarsky et al.l 
(2006c) showed that the best observational targets are the dwarf 
spheroidals (Ursa Minor, Draco, etc). Indeed, these objects are 
X-ray quiet, while at the same time one expects the DM decay 
signal from them, comparable to wh at comes from galaxy cl us- 
ters. Because at the time of writing of lBovarskv et al.l d2006d) no 
public data were available for these dwarf spheroidals, the ob- 
servations of the core of Large Magellanic Cloud (LMC) were 
used to produce the strongest rest rictions on parameters o f the 
sterile neutrino. It was stressed in Boyarsky et al. (2006c) that 
other dwarf satellite galaxies should be studied as well, in or- 
der to minimize uncertainties related to the DM modeling in 
each single object. In this paper we continue studies of the 
dwarf satellites of the MW by analyzing the data from XMM- 
N ewton observation of Ur sa Minor and confirm the restrictions 

of lBovarskv et al l (l2006d) 

It was also shown in IBovarskv et al.l (l2006d) that the im- 
provement of the results from MW DM halo can be achieved 
by using longer exposure data (notably, longer exposure of the 
closed filter observations). In this paper, we improve our restric- 
tions, coming from the MW DM halo by using the bla nk sky 
dataset with better statistics from Nevalaine n et al.l (120051) . 



2. DM with radiative decay channel 

Although throughout this paper we are talking mostly about the 
sterile neutrino DM, the results can be applied to any DM par- 
ticle that possesses the monoenergetic radiative decay channel, 
emits photon of energy E y and has a decay width V. In the case 
of the sterile neutrino (with mass below that of an electron), 
the r adiative decay channel is into a photon and active neu- 
trino iPal & Wolfensteinl (1982). As the mass of an active neu- 
trino is much lower than keV, E y = in this case. The width 
r of radiative deca y can be expressed (IPal & Wolfensteinlll982t 
Barg er et"ai1l 19951) in terms of mass M s and mixing angle 6 via 



= 9aG F sin 2 26M 5 „ i.38 X io- 22 sin 2 (26)) 
1024tt 4 



M, 



keV 



■(1) 



(The notation sin 2 (2#) is used traditionally, although in all re- 
alistic cases 8 « 1), The flux of the DM decay from a given 
direction is given by 

PdmOO 



e 7 r 



fov cone 



An\D L + if 



T d 3 r 
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Here Dl is the luminous distance between the observer and the 
center of the observed object, PdmW is the DM density, and 
the integration is over the DM distribution inside the (truncated) 
cone - solid angle, spanned by the field of view (FoV) of the 
X-ray satellite. If the observed object is farQ then Eq. (f2| can be 
simplified: 



DM 



DM ^7 

4nD 2 L M s 



(3) 



2 Namely, if luminosity distance D L is much greater than the charac- 
teristic scale of the DM distribution pou(r). 



where M^ M is the mass of DM within a telescope's field of view 
(FoV). Equation (0 can be rewritten again as 
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In the absence of a clearly detectable line, one can put an upper 
limit on the flux of DM from the astrophysical data, which will 
lead via Eq. (01 to the restrictions of parameters of the sterile 
neutrino M s and 8. 



3. Restrictions from the blank sky observation 

3.1. Modeling the DM halo of the MW 

As shown in the previous section, one needs to know the distri- 
bution of the DM to obtain the restrictions on parameters of the 
sterile neutrino. In the case of nearby objects (including our own 
Galaxy and dwarf satellites from the local halo), the DM dis- 
tribution can be deduced e.g. by using the rotati on curves of the 
stars in the galaxy. Here we follow the analysis of Boyarsky et al.l 
(2006c). Various DM profiles, used to fit observed velocity dis- 
tributions, differ the most in the cen ter of a distribution. I n 
the case of the MW we choose, as in IBovarskv et al.l (|2006c), 
to use the observations away from the center, to m i nimiz e 
this uncertainty. In p articular, in Refs. iKlvpin et al.l d2002l) : 
iBattaglia etaf] (I2005I) it was shown that the DM halo of the 
M W can be described by the Navarro-Frenk- White (NFW) pro- 



file dNavarro et 



.cribe d by 
ai1 ll997l) 



PnfwW 



r(r + r s ) 2 



(5) 



with parameters, given in Table Q~Q The relation between virial 
parameters and p s , r s are given in the Appendix lAl Quoted halo 
parameters provide DM decay flux (from the directions with 
<p > 90°) consistent within ~ 5% wit h the one, given by Eqs. (O- 
©. Only "maximal disk" models in Klvpi n~et al.l d2002l) would 
provide 30 - 50% weaker restrictions; howeve r, these models 
are highly implausible, see Klyp in et al.l ( 20021). Similarly, tak - 
ing the lower limit for the virial mass of Battagl ia et al.l y005), 
one would obtain 25% weaker restrictions than the ones, pre- 
sented in this paperQ 

To compare the results from different (e.g. cuspy and cored) 
profiles, we can also describe the DM distribution in the MW via 
an isothermal profile: 



PisoM 



1 



4-nG N r 2 + r 2 



(6) 



The DM flux from a given direction <p into the solid angle £2f ov <k 
1, measured by an observer on Earth (distance r 8 kpc from 
the galactic center), is given by 



dmv 



Lq f§ + arctan(^f^) , cos > 



1 = — X 



R 1 arctan ( R ,, ) , cos0<O 



(7) 



3 According tolKlvpin et al. (2002), the choice of e.g. the Moore pro- 
file lGhigna e t al. (2000) or a generalization thereof, as compared to the 
NFW profile, would change the results by < 1% for r < 3kpc. As we are 
using observations away from the center, this difference is completely 
negligible, so we choose to use the NFW profile. 

4 When quoting results of IKlvpin et al.l d2002l) . we do not take the 
effects of baryon compression on DM into account. While these ef- 
fects make DM distribution in the core of the MW denser, they are hard 
to compute precisely. Thus the values we adopt give us a conservative 
lower bound on the estimated DM signal. 
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Table 1. Best-fit parameters of NFW model of the MW DM halo. 
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References 


Myir [M Q ] 


fvir [kpc] 


Concentration 


r s [kpc] 


p s [M Q /kpc 3 ] 


Klvpin et al. (2002), favored models (Aj or B\) 


1.0 x 10 12 


258 


12 


21.5 


4.9 x lO" 


Battaglia et al. (2005) 


0.8+, 1 , 2 x 10 12 


255 


18 


14.2 


11.2 x 10 6 



Here Lq = ^^"cn anc ' ^ = + r o sm2 ^' Angle ^ s reia ted to 
the galactic coordinates (b, I) via 

cos <p = cos b cos / . (8) 

Thus, the galactic center corresponds to <p — 0°, and the anti- 
center (p = 180° and the direction perpendicular to the galactic 
plane t o = 90°. 

In iBovarskv et alj d2006d) . the following parameters of 
isothermal profile were chosen: v/, = 170 km/sec and r c = 4 kpc. 
One can easily check (using Table Q] and Eqs. dA.5t - dA.6b in 
Appendix [Ab that, in the directions (p > 90°, the difference in 
predicted DM fluxes between the NFW model with parame- 
ters, given in Table[TJand isothermal model with parameters just 
quoted are completely negligible (less than 5%)Q 

3.2. XMM-Newton PN blank sky data 

To examine the Milky Way halo, we used the double-filtered, 
single+double event X MM-Newton PN blank sky data from 
iNevalainen et al.l ((2005), whi ch is a collection of 18 blank sky 
observations (see Table 2 in INevalainen et_ al.l d2005l) for their 
observation IDs, positions, and exposures) The exposure time 
of the co-added observations is 5 47 ks. We used a combi nation 
of closed-filter observations from INevalainen et al.l (|2005) (total 
exposure time 145 ks) to model the background of XMM-Newton 
PN instrument separately. The data has been filtered using SAS 
expression "flag==0", which rejects the data from bad pixels and 
CCD gap regions. After removing the brightest point sources, the 
total accumulation area is 603 arcmin 2 . 

Based on the >10 keV band count rates of the blank sky 
and the closed filter data, we normalized the closed-filter spec- 
trum by a factor of 1.07 before subtracting it from the blank 
sky spectrum. The remaining sky-background spectrum consists 
mainly of the Galactic emission and the cosmic X-ray back- 
ground (CXB) due to unresolved extragalactic point sources. 
We modeled the Galactic emission by a non-absorbed MEKAL 
model with Solar abundances. For the CXB emission, we used 
a power-law model modified at the lowest energies by Galactic 

5 Ref . lAbazaiian & Koushiappas (2006) claim that the MW results 
of Refs. Bovars kv et al.l d20 06c): Rie mer-S0rensen et al.l d2006T) are un- 
certain by about a factor of 3. This conclusion was based on the range 
of viria l masses of t h e MW DM halo M vir = (0.7 - 2.0) x 10 12 M 
in Ref. iKlypin et alj d2002l). H owever, as just demonstrated, the au- 
thors of Bovarskv et all d2006c) have chosen parameters of DM halo 
conservatively. The flux they used , corresponded to the favored mod- 
els A\ or Bi in lKlvpin et aild2002h . with M vir ~ 1.0 x 1O 12 M . These 
models provide the lowest bound on the derived flux of DM decay 
(if one does not take into acc ount the highly im plausible "maximum 
disk" (A 2 or fi 2 ) models of IKlypin et all d2002h ). Even in the lat- 
ter c ase, the DM flux will b e 30 - 50% lower than the one used in 
work IBovarskv etldl d2006ch. Th erefore, parameters of the MW DM 
halo from lBovarskv et all {2006c) provide the conservative estimate so 
we use them in our work as well. 

6 We processed the blank sky data with newer SAS distribution, 
xmmsas_28050815_1803-6. 5.8, and obtained slightly different expo- 
sure times than those in the public data. 




0.6 - 



1.0 2.0 5.0 

Energy [keV] 

Fig. 1. Upper panel: The blank sky data of INevalainen et aTl 
d2005l) after subtraction of the closed filter data. The black 
crosses show the data points, used for modeling the Galactic 
emission and CXB. The blue crosses show the data points ex- 
cluded from the fit. The best-fit model is shown with the red line. 
The plotted spectrum is binned by a minimum of 20000 counts 
per channel, different from that used in the actual fit (see text). 
The error bars include the 5% systematic uncertainty used in the 
analysis. Lower panel: The ratio of the data-to-model values in 
each channel used in the fit. 



absorption with the value of Nh fixed to its exposure-weighted 
average over all blank sky observations (Nh = 1.3 X 10 20 cm -2 ). 

The variable Galactic emission and geocoronal Solar wind- 
charge exchange emission (see e.g. IWargelin et alj|20 04) com- 
plicate the modeling at the lowest energies. The remaining cali- 
bration uncertainties further complicate the analysis in the low- 
est energies (see e.g. INevalainen et ai1l2006l) . Thus, we omitted 
the channels below 0.8 keV. At energies above 7 keV, the particle 
background dominates and the total flux is very sensitive to the 
background normalization. We thus excluded channels above 7 
keV. 

The data are not well-described in the 1.45-1.55 and 5.8- 
6.3 keV bands with the above model. These deviations proba- 
bly originate from the variability of the instrumental Al and Fe 
line emission. In order to minimize the effect of the instrumen- 
tal problems, we excluded these bands when finding the best-fit 
sky background model (see below). Also, to account for possi- 
bly remaining calibration inaccuracies, we added a systematic 
uncertainty of 5% of the model value in each bin in quadrature 
to the statistical uncertainties. 

We binned the spectrum using a bin size of 1/3 of the energy 
resolution and fitted the data using models and channels as de- 
scribed above. The best-fit (reduced^- 2 = 1.03 for 153 degrees 
of freedom) model agrees with the data within the uncertainties 
(see Fig. [TJ, yielding a photon index of 1.50+0.02 at 1 <x confi- 
dence lev el (see Fig[Tjl, consistent with similar analyses based on 
Chandra dHickox & Markevitch[2006l) and XMM-Newton MOS 
instrument dDe Luca & Molendil20 04). The best-fit temperature 
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Fig. 2. Exclusion plot based on the blank sky observations. Bins 
2.9 - 3.1 keV and 11.6 - 12.6 keV are excluded. 



of the MEKAL component used to mode l the Galactic emission 
is 0.1 9±0.01 keV, consistent with e.g. iHickox & Mark evrtohl 
(120061) . 

We then evaluated the level of possible DM flux above the 
background model allowed by the statistical and systematic un- 
certainties in each channel. For this, we modified the above best- 
fit model by adding a narrow (width = 1 eV) Gaussian line to it. 
We then re-fitted the data, fixing the Gaussian centroid for each 
fit to the central energy of a different channel. In these fits we 
fixed the above continuum model parameters to the best-fit val- 
ues and thus the Gaussian normalization parameter is the only 
free parameter. We used the fits to find the upper 3<x uncertainty 
of the Gaussian normalization, i.e. the allowed DM flux. Note 
that here we included the channels 1.45-1.55 and 5.8-6.3 keV 
(excluded above when defining the sky background model). The 
background is oversubtracted in the channels at 1.45-1.55 keV 
and ~6.0 keV (see Fig. [D, which would formally require neg- 
ative normalization for the Gaussian. However, we forced the 
normalization to be positive and thus obtained conservative up- 
per limits in these energies. 

Finally, we converted the upper bound obtained for the 
flux per energy bin to the restrictions on M s and sin 2 29, us- 
ing Eqs. |[l), © (we use exposure weighted average of DM 
fluxes (Eq. |7J from all the observations, constituting the blank 
sky dataset). This corresponds to the average "column density" 
1.22 x 10 28 keV/cm 2 . The results are shown in Fig. [2] 

At energies above E = 5 keV, the instrume ntal background 
of PN dominates over the sky background (c.f. lNevalainen et all 
120051) . Therefore, the accuracy of the co-added closed filter spec- 
trum in predicting the particle background in the blank sky ob- 
servations becomes essential. We estimate this accuracy using 
the variability of the individual close d-filter spectra in the 0.8- 
7.0 keV band (Nev alainen et a l. 2005) and propagate it by vary- 
ing the normalization of closed filter data by +5% and repeating 
the above analysis. This leads to a factor of 3 change in the re- 
sults at M s ~ 14 keV (see FIG. [3). Therefore, for E > 5 keV we 
choose the more conservative normalization (see FIG.|7]below). 



4. Restrictions from observations of Ursa Minor 



It was argued in [Boyarsky et alj d2006d) that dwarf satellite 
galaxies should provide the best restrictions, based on their high 
concentration of DM and low X-ray signal. At the moment 
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Fig. 3. Dependence of the results on closed filter normalization. 
Red (solid) and green (long-dashed) lines are the same as in 
Fig. [2] 



of writing of Boyarsky et al.l d200 6c). no public data on pre- 
ferred dwarf satellites were available, therefore the observation 
of core of LMC were used as a demonstration. Recently, the 
Ursa Minor dwarf (UMi) was observed with XMM-Newton (obs 
IDs.: 0301690201, 0301690301, 0301690401, 0301690501, ob- 
served in August-September 2005)Q Unfortunately, most of 
these observations are strongly contaminated by background 
flares and the observations have very small exposure times. 
Below we present the analysis of only one observation (obsID: 
0301690401), which "suffered" the least from background con- 
tamination. 



4.1. DM modeling for UMi 

The DM distribution in UMi has a co red profile (see e.g. 
i Klevna et ai]|2003l IWUkinson et al.ll2006t iGilmore et al.ll2006l 
2007)0 We adopt the following parameters of isothermal pro- 
file (|6]l_jOT_lJMj^V£_ = 23 km/sec, r c = 0.1 kpc (see 
e.g. lWilkinson et al.i 2006)0 We adopt the distance to UMi D L = 



7 We are very grateful to Prof. T. Maccarone for sharing this data with 
us before it became publicly available through the XMM data archive. 

8 As discussed in Section [3"Tl the estimates for DM flux do not vary 
significantly if one uses NFW instead of the isothermal DM density 
profile. In the case of UMi, the cored (isothermal) profile will clearly 
produce a more conservative estimate than will the cuspy NFW pro- 
file. In deed, taking NFW parameters for UMi from the recent paper lWu 
( 2007) gives a ~ 20% higher estimate for the DM mass within the FoV. 

9 For the detailed stud ies of mass distribution in dwarf spheroidals, 
see lGilmore et al.1 d2007h . We are grateful to Prof. G. Gilmore for shar- 
ing the numbers with us before their paper became available. The sta- 
tistical uncertainty in determining these numbers is below 10%. The 
systematic uncertainties are much harder to estimate. One of the ma- 
jor sources of the systematic errors comes from violation of the main 
assumptions of the method: deviation from equilibrium and from the 
spherical distribution of matter in a galaxy. In other known examples 
it provides a factor of 2 uncertainty, which should be a conservative 
estimate in the case of UMi, as it is rather spherical. Another typical 
uncertainty - determination of the mass of the stars - is not important 
for UMi, as it has a very high mass-to-light ratio. 
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2 4 6 8 

Time from PN obs. stort [ks] 

Fig. 4. Hard band light curve for the UMi observation 
0301690401. The crosses show the PN >10 keV band rate of 
the full FOV in Iks bins. The dashed line shows the average, 
when excluding first 2 time bins. The dotted line shows the 
corresponding quiescen t value in the co-added blank sky data 
dNevalainen et al]|2005l) . 



66kpc dMateoll998l) . The DM mass within the circular FoV with 
the radius rf 0V , centered at the center of the galaxy is given by 



ZLtm 



(9) 



In our case, the radius of FoV is 13.9', which corresponds to 
rfcw = 0.27 kpc (i.e. about 3r c ). Therefore 



Ml 



i'ov 



3.3 x 1O 7 M for r fov = 0.27 kpc. 



(10) 



Using Eqs. ( fTOb and (0), one can compute the expected DM flux 
from UMi: 



DM 



: 4.79 



keV 



cm' ■ sec 



keV 



sin 2 (20). 



(ID 



4.2. PN data analysis 



We processed the Ursa Minor observation 0301690401 us- 
ing epchain version 8.56 and filtered the event file with 
SAS expressions "PATTERN<=4" and "FLAG==Q". We applied 
the blank sky-b a sed X MM-Newton background method of 
iNevalainen et al.l d2005l) for Ursa Minor. The > 10 keV band 
light curve from the full FOV (FIG. |4]l shows that the count rate 
in observation 0301690401 (excluding first 2ks) exceeds that of 
the blank sky quiescent average by 25%. This level is higher, 
but close to what is used in the blank sky accumulation (+ 20% 
filtering around the quiescent level). Thus we accepted the data 
from all instants after the initial 2ks, and we approximated th e 
background uncertainties with those in Nevalain en et al.l (|2005). 

The hydrogen column density in the direction of Ursa Minor 
is small (~Nh = 2 x 10 20 crrr 2 ) and consistent with the variation 
in the blank sky sample. Thus we can also apply the blank sky 
background method to channels below 2 keV. 

As noted in the above XMM-Newton blank sky study, the > 
10 keV band-based scaling of the background only works up to a 
factor of 1 . 1 , beyond which the background prediction becomes 
worse. Furthermore, the correlation of background rates in the > 
10 keV band is very poor with the rates below 2 keV band. Thus, 
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Fig. 5. Flux from UMi (obs. ID 0301690401). Energy bins have 
the width of twice the spectral resolution. Shown are the 1, 2, 
and 3<x errors. One can see that, above 2 keV, flux in most en- 
ergy bins is zero within Icr limits (blue crosses) and for the rest it 
is zero within 2<x limits (green crosses). Similarly, below 2 keV 
black, cyan and yellow crosses represent 1,2, and 3<x error corre- 
spondingly. The solid black line represents the 3<x upper bound 
on total flux in a given energy bin, which we use to put the limit 
on DM parameters. 



in order to achieve the best possible background prediction ac- 
curacy, we scaled the blank sky background spectrum by a factor 
of 1.1 at channels above 2 keV, and we applied no background 
scalingat lower energies. 

We removed this scaled background spectru m from the Ursa 
Minor spectrum (see FIG. |5}. As shown in INevalainen et al.l 
(2005), the background accuracy is worse at lower energies. We 
used those estimates to propagate the background uncertainties 
at 1 cr confidence level to our results by examining how the 
results change when varying the 0.8-2.0 keV and 2.0-7.0 keV 
band background by 15% and 10%, respectively. 

4.3. Ursa Minor data and restrictions on the sterile neutrino 
parameters 

The X-ray spectrum of UMi is similar to that of LMC: above 
2 keV the flux is zero within statistical limits (see FIG. [5]). (Of 
course, the data set has rather low statistics: after the cleaning 
of flares the UMi observation only contains 7 ks). Therefore, 
for such data, we utilized the "total flux" method. Namely, we 
restricted the DM flux in the given energy bin to be bounded 
from above by the measured total flux in this energy bin plus its 
3<x uncertainty (FIGHJ. As each energy bin has a width of twice 
the FWHM at a given energy, the flux from a DM line would not 
"spill" into nearby bins. Using Eq. (11 It . we find the restrictions 
on the sterile neutrino parameters, represented in FIG. [6] 

These restrictions should be compared with those, ob- 
tained from another satelli te galaxy - Large Magellanic 
Cloud dBovarskv et al.ll2006d) . As one clearly sees from FIG. [6] 
in spite of the low exposure time, it is fully consistent with 
th e earlier bound s from LMC, thus confirming the results 
of iBovarskv et al.l d2006d) . Improvement of the exposure for 
UMi observations should, presumably lead to the improvement 
of results (at least for energies above E > 2 keV). For example, 
for a 100 ksec observation, we expect the results to improve by 
roug hly a factor a/100 ks/7 ks * 3.77. 
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Fig. 6. Exclusion from UMi (red solid line), as compared to 
LMC (green long-dashed line). The blue short-dashed smooth 
curve shows hypothetical restrictions from UMi observations 
with 100 ksec exposure. 
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Fig. 7. Results and comparison with previous bounds (a region 
of parameter space above a curve is excluded) 



5. Results 

5.1. Restrictions from the blank sky data 

By analyzing the blank sky data set with better stati s tics, w e 
improved on the prev i ous results of Boyarsk y et al.l (|2006c), 
lRiemer-S0rensen et all (l2006l) . IWatson et al.l(l2006h (by as much 
as the factor of 10 for M s » 3.5 keV and by the negligible 
amount for M s > 11 keV). The result is shown in Fig. [7] in 
red solid li ne. The best previous b ounds are also shown: bound 
from LMC dBoyarsky et al. 2006c) wi th a bl ue short-dashed line 
and bound from M31 (Wats on et al.l |2006h with a dotted ma- 
genta line. We see that in the region 3.5 keV < M s < 11 keV 
the new blank sky data improves on previous results. These re- 
sults can be converted (using Eq. ([T])) into restrictions on the 
decay rate F of any DM particle, that possesses radiative decay 
channel and emits a photon E y (see FIG. [8]). Our results provide 
more than an order of magnitude improvement ov er similar re- 
strictions derived in Riemer-S0re nsen et alj d2006l) (which used 
the Chandra blank sky backgro und), as one can clearly see b y 
comparing FIG. [8] with FIG. 2 in Rie mer-S0rensen et al.l ((2006), 
where the exclusion plot is above F = 10~ 26 sec -1 line for all 
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Fig. 8. Restrictions on parameters of any DM particle with the 



radiative decay width F, 
region is excluded. 



emitting photon of energy E y . Shaded 



energies. (In lRiemer-S0rensen et alj (|2006) the restriction were 
made, based on the total flux of Chandra satellite, without sub- 
traction of the instrumental background, which explains a much 
weaker restrictions). 

The empirical fit to the MW data is given by the following 
expression: 



sin 2 (20) < 2.15 x 10~ 7 (— ) 
VkeV/ 



-3.45 



(12) 



5.2. Restrictions from Ursa Minor dwarf 



Restrictions from XMM observation of UMi are shown in FIG. [7] 
by the green long -dashed line. These results are slightly weaker 
than LMC or M3 1 results, which is due to the very low statistics 
of the UMi observation. Improvement of the statistics should 
lead to improvement of the current bound (as shown on the 
FIG. |6]l. These results confirm the recent claims (Boyarsky et al., 
2006c) that dwarfs of the local halo are promising candidates for 
the DM-decay line search and, as such, should be studied dedi- 
catedly. 

In searching for the DM signatures, it is important to under- 
stand that the uncertainties of the DM modeling for any given 
object can be large, and therefore it is important to study many 
objects of given type, as well as many different types of objects 
(where DM distributions are deduced by independent methods). 
To this end, although UMi data does not provide any improve- 
ment over existing bounds, it makes those bounds more robust 
as the existence of DM in UMi is deduced by independent obser- 
vations, and the rotation curves of UMi are measured quite well, 
since it has less perturbed dynamics, compared to e.g. LMC. 

6. Discussion 

In this paper we continued to search for the best astrophysical 
objects, from the point of view of restricting parameters of DM 
particles with the radiative decay channel. Several comments are 
in order here. 

(1) Although throughout this paper we have spoken about the 
sterile neutrinos and restricted their parameters (namely, 
mass, and mixing angle), the constraints can be readily con- 
verted into any other DM candidate that possesses a radia- 
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(2) 



tive decay channel. [3 Then the restrictions are formulated 
on the decay rate T as a function of energy of the emitted 
photon E y . The results then can be presented in the form of 
an exclusion plot, presented in FIG. [8] 
Clearly, if one could relate parameters of the sterile 
neutrino with their relic abundance Q s , this would 
allow one to put an upper limit on the mass of the 
sterile neutrino. Unfortunately, s uch a computation 
is str o ngly model-dependent. In [D odelson & Widrow 
( 1994), Dolgov & Hansen! d2002l) . lAbazaiian et alJ d2001al) . 
lAbaza iian (20061). the relic abundance of the sterile 
neutrinos was computed in a simple model with only 
one sterile neutrino, assuming the absence of the sterile 
neutrinos above the temperatures ~ 1 GeV. Yet, even 
the computation in this simplest model is subject to a 
number of uncertaintie s dShi & Fullerlll999HBovarskv et"ai1 
2006all3: lAsaka et al.1 l2006allbt IS haposhni kov & Tkachevl 
20061). In partic ula r, in iDodelson & Widrowl (|1994 | ) , 
Dolgov & Hansen! (|2002|) . lAbazaiian et al.1 ( 1200 h 



Abazajian (2006), two assumptions were made: (i) the 



absence of heavy particles, whose decay can dilute the relic 
abundance, (ii) the absence of lepton asymmetries. In addi- 
tion, simplifying assumptions about dynamics of hadrons at 
tempe ratures (9(150) MeV were used. Recently, Asak a et ail 
(2007) performed this computation from the first principles, 
showing that the uncertainty due to QCD effects (between 
minimal and maximal values of sin 2 (2#) for given M s ) is 
about a factor of 8. 

Taking away the assumptions about the absence of the sterile 
neutrinos above the temperatures ~ 1 GeV makes any mix- 
ing angle possible. For example, the DM neutrinos can be 
created due to the inflaton decay Shap oshnikov & Tkachevl 
(2006). Therefore, in this work we chose not to derive an 
upper bound on the mass of the sterile neutrino. 
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Appendix A: Determining parameters of NFW 
profile 

Using the data on rotation curves, one usuall y obtains the follow - 
ing parameters of DM distribution (see e.g. Kly pin et al.l d2002)): 
virial mass M wu , virial radius r v i r , and concentration parameter 
C. They have the following relation with the parameters of NFW 
profile © r s and p s \ 



M v 



_ ' vir _ 

where in terms of function fix) 
/(x)=log(l+x)- 



1 + X 



(A.l) 



(A.2) 



10 For earlier works, discussing c osmological and astro p hysica l 
effects of decaying DM, see e.g. Ide Ruiula & Glashowl ( l980); 
IBerezhiani et alJ il987h: l5 oroshkevic h et alJ dl989h : iBerezhiani et alJ 
( 1990); Berezhiani & Khlopov ( 1990) . The extens ive review of the re- 
sults can also be found in the book Khlopov] jl997l) . 



one obtains the mass within the radius r. 

f(r/r s ) 



M(r) = M V] 



f(C) 



(A3) 



If DM distribu tion in the Milky Way is described by the NFW 
model (as in lBattaglia et al.l 120051 iKlypin et al.ll2002l) . the flux 
from a direction (f> is given by 



77 NFW / i \ 



f'OY 



8;r 



J"t fepNFw( -\J> 



ri+z 2 + 2r Q z cos ( 



(A.4) 



(notations are the same as in Eqs.(|7|i-([8|l). Let us consider 
two cases, when the integral in ( 1A.41 i can be easily computed. 
Namely, we have for <p = 180° 



7 NFW 



(180°) = 



fov 



8;r 



-p s r s 



log(l + J.) £_ 



(A.5) 



and for <f> = 90° 



7 NFW 



(90°) 



FQ 



fov 



8;r 



-p s r s 



■1 -los 



2r s 



- 2 (3 



log 



2r s 



+0 



(in the latter case the analytic expression is too complicated, so 
we present Taylor expansion for the case r <K r s ). 

References 

Abazajian, K. 2006, Phys. Rev., D73, 063506 

Abazajian, K., Fuller, G. M., & Patel, M. 2001a, Phys. Rev., D64, 023501 
Abazajian, K., Fuller, G. M., & Tucker, W. H. 2001b, Astrophys. J., 562, 593 
Abazajian, K. & Koushiappas, S. M. 2006, Phys. Rev., D74, 023527 
Aguilar, A. et al. 2001, Phys. Rev., D64, 112007 

Asaka, T, Blanchet, S., & Shaposhnikov, M. 2005, Phys. Lett., B631, 151 
Asaka, T, Kusenko, A., & Shaposhnikov, M. 2006a, Phys. Lett., B638, 401 
Asaka, T, Laine, M., & Shaposhnikov, M. 2006b, JHEP, 06, 053 
Asaka, T, Laine, M., & Shaposhnikov, M. 2007, JHEP, 01, 091. 
Asaka, T. & Shaposhnikov, M. 2005, Phys. Lett., B620, 17 
Barger, V. D., Phillips, R. J. N, & Sarkar, S. 1995, Phys. Lett., B352, 365 
Battaglia, G. et al. 2005, Mon. Not. Roy. Astron. Soc, 364, 433 
Berezhiani, Z. G. & Khlopov, M. Y. 1990, Sov. J. Nucl. Phys., 52, 60 
Berezhiani, Z. G, Vysotsky, M. I., & Khlopov, M. Y. 1987, Sov. J. Nucl. Phys., 
45, 1065 

Berezhiani, Z. G., Vysotsky, M. I., Yurov, V. P., Doroshkevich, A. G, & Khlopov, 

M. Y. 1990, Sov. J. Nucl. Phys., 51, 1020 
Biermann, P. L. & Kusenko, A. 2006, Phys. Rev. Lett., 96, 091301 
Bode, P., Ostriker, J. P., & Turok, N. 2001, Astrophys. J., 556, 93 
Boyarsky, A., Neronov, A., Ruchayskiy, O., & Shaposhnikov, M. 2006a, 

MNRAS, 370, 213 

Boyarsky, A., Neronov, A., Ruchayskiy, O., & Shaposhnikov, M. 2006b, 

Phys. Rev. D, 74, 103506 
Boyarsky, A., Neronov, A., Ruchayskiy, O., Shaposhnikov, M., & Tkachev, I. 

2006c, Phys. Rev. Lett., 97, 261302 
de Gouvea, A. 2005, Phys. Rev., D72, 033005 
De Luca, A. & Molendi, S. 2004, Astron. Astrophys., 419, 837 
de Rujula, A. & Glashow, S. L. 1980, Phys. Rev. Lett., 45, 942 
Dodelson, S. & Widrow, L. M. 1994, Phys. Rev. Lett., 72, 17 
Dolgov, A. D. & Hansen, S. H. 2002, Astropart. Phys., 16, 339 
Doroshkevich, A. G, Khlopov, M. I., & Klypin, A. A. 1989, MNRAS, 239, 923 
Ghigna, S., Moore, B., Govemato, E, et al. 2000, ApJ, 544, 616 
Gilmore, G, Wilkinson, M., Kleyna, J., et al. 2006, astro-ph/0608528. 
Gilmore, G, et al., 2007, astro-ph/0703308. 

Goerdt, T, Moore, B., Read, J. I., Stadel, J., & Zemp, M. 2006, MNRAS, 368, 
1073 

Hansen, S. H., Lesgourgues, J., Pastor, S., & Silk, J. 2002, MNRAS, 333, 544 
Hickox, R. C. & Markevitch, M. 2006, ApJ, 645, 95 
Hidaka, J. & Fuller, G. M. 2006, astro-ph/0609425. 

Khlopov, M. Y. 1997, Cosmoparticle Physics (World Scientific Pub Co Inc) 
Kleyna, J. T, Wilkinson, M. I., Gilmore, G., & Evans, N. W. 2003, ApJ, 588, 
L21 

Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 
Kusenko, A. 2006a, astro-ph/0609375. 



8 



A. Boyarsky et al.: Dark matter constraints from the Galaxy and Ursa Minor 



Kusenko, A. 2006b, Phys. Rev. Lett., 97, 241301 
Mateo, M. L. 1998, ARA&A, 36, 435 

Navarro, J. R, Frenk, C. S., & White, S. D. M. 1997, Astrophys. J., 490, 493 
Nevalainen, J., Bonamente, M., & Kaastra, J. 2006, astro-ph/06 10461. 
Nevalainen, J., Markevitch, M., & Lumb, D. 2005, Astrophys. J., 629, 172 
Pal, P. B. & Wolfenstein, L. 1982, Phys. Rev., D25, 766 
Riemer-S0rensen, S., Hansen, S. H., & Pedersen, K. 2006, ApJ, 644, L33 
Riemer-SOrensen, S., Pedersen, K., Hansen, S. H., & Dahle, H. 2006, astro- 
ph/0610034. 

Seljak, U., Makarov, A., McDonald, P., & Trac, H. 2006, Phys. Rev. Lett., 97, 
191303. 

Shaposhnikov, M. & Tkachev, I. 2006, Phys. Lett., B639, 414 
Shi, X.-d. & Fuller, G. M. 1999, Phys. Rev. Lett., 82, 2832 
Stasielak, J., Biermann, P. L., & Kusenko, A. 2006, astro-ph/0606435. 
Strigari, L. E. et al. 2006, astro-ph/0603775. 
Strumia, A. & Vissani, F., hep-ph/0606054. 

Viel, M., Lesgourgues, J., Haehnelt, M. G., Matarrese, S., & Riotto, A. 2005, 

Phys. Rev., D71, 063534 
Viel, M., Lesgourgues, J., Haehnelt, M. G, Matarrese, S., & Riotto, A. 2006, 

Phys. Rev. Lett., 97, 071301 
Wargelin, B. J., Markevitch, M., Juda, M., et al. 2004, ApJ, 607, 596 
Watson, C. R., Beacom, J. R, Yuksel, H, & Walker, T. P. 2006, Phys. Rev., D74, 

033009 

Wilkinson, M. I. et al., 2006, astro-ph/0602186 
Wu, X., 2007, astro-ph/0702233 



